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ABSTRACT 

The present work asseses the reliability and accuracy of such a 
formulation for bending analysis of laminated plates under different physical 
conditions and using various material parameters. The FEM adaption of this 
formulation has been carried out by using , isoparametric, Lagrangian quadratic 
elements . Transverse shear stresses have been evaluated using the constitutive 
relations as well as the equilibrium equations . A large number of experiments 
were carried out to establish the reliability of this theoretical formulation and 
its FEM adaption. The results have been compared with those from the First Order 
Shear Defcrmation Theory and CPT in order to determine the domain of their 
applicability. A FORTRAN 77 program has been developed for the computer 
implementation of this formulation . 



CHAPTER 1 


INTRODUCnON AND LITERATURE SURVEY 


l.i fritroducticn : 

Now a days fiber reinforced composite materials are used extensively 
as structural members in high to low techwTology areas. Their preferred choice is 
due to their superior mechanical properties. The numerous combinations of 
ccnstituent materials and numerous options in fiber cx'ientation and stacking 
sequence permit a wide scope in composite structural design. A composite material 
offers a considerable weight reduction due to its high strength to weight and 
stiffness to weight ratio. Usually other properties such as fatigue life, corrosion 
and wea^ resistance, temperature dependent behavicw' etc. are improved 
considerably. 

Fibrous composite material is a two phase system consisting of high 
strength fiber as main load bearing member, embedded in a matrix. Some of the 
mcff'e conBTion type of fibers that are in use are graphite, glass, boron, silicon 
carbide, beryllitm and carbon. Epoxy phenols, polyester, altminum, titanium and 
epoxies are a few common matrix materials being used at present. When these two 
(^ase systan are sL±ijected to various loading conditions, the matrix transfer 
stresses to the high strength fibers through inter facial bonds between fiber and 
matrix. In addition to this matrix also helps in keeping fibers in their appropriate 
positicn arad protects them from environmental effects and provides flexibility to 
Fiba' in addition to reinforcing matrix, tends to arrest any crack 
development rrarmal to its direction . For the sake of mathematical modeling a 
lamina is caisicte’^ to be a homogeneous, 2— dimensional, orttotropic layer . In 
general a laminate i«#»ich is a ccmbination of several laminae placed over each 
otha' is ccwisidered heterogeneous tlr'ou^ the thickness and generally 
anisotrc«)ic. 

Delamination has become a problem of great concern. Initiation and 
»’x>wlh of this kind of failure mode is ckie to inter— laminar stresses acting 
throjoh the thickness. Strength requirements are governed by in— plane stresses. 
TTkis a tlreory i^hit±i can predict all this stresses accxiralely becomes irmoerative 
fo" efficiait and reliable desig^i and analysis of structural mentier made up of 
ccxnposite materials. 



Due to the hi#i difference in properties between matrix material and 
fiber material, there is a greater difference in in— plane Young's modulus and 
transverse shear modulus as compared to isotropic materials. Hence classical 
lamination theory which is an extension of classical plate theory can no longer 
pr«iict ccM^rect responses in the case of multi layered composite structure. 
Ftff'ther investigations in this regard lead to first order shear deformation theory 
i4iich assumes constant shear rotation and requires use of shear correction 
factor. The prediction of those shear correction factor in the case of anisotropic 
material is c^ite involved and is problem dependent. Also this theory does not take 
into accotnt the warping of transverse cross section which takes place specially 
in a thick and weak core of a sandwich construction. Further the theory neglects 
the effects of transverse normal stress / strain. These limitations of first order 
shear deformation theory prompted for the development of a more refined theory 
t#)ich considers realistic parabolic variation of transverse shear stresses 
through the plate thickness, warping of the transverse cross section and the 
complete Hook's law. 

12 Literature Review : 

The literature dealing with flexural analysis of laminated composite 
plates can be categorized among three groups. The first being called as classical 
lamination (thin plate ) theory i4iich is an extension of classical plate theory to 
laminated conn^osite plates. This theory ignores the effects of shear straind, 
normal strains, and normal stresses in transverse direction and hence it fails to 
predict responses acccrately in the case of laminated plates. These limitatiwis of 
(XT have lead to the development of first order shear deformation theory. The 
literatir«s relating this first order shear deformation theory comprises the 
second group. The literatures dealing with the most refined theories based on 
either 3— dimansional ac«)roach or 2 —dimensional approach with higher order 
displacement models i4iic±> take accoint of realistic parabolic variation of 
transverse shear stress throu^ the plate Wrickness ( thus reoJiring ra use of 
shear cor^stion factor ) and transverse cross section warping, ccMi^oises tl« 
Wiird groL«3. 

The first cowlete laminated anisotrcpic plate theoHj is attributed to 
Reiasnor and Stavs^y^. It inoludas oo^nling between in— plana stretching and 



tranffVQrsfii bandingf for unsymmatrioally laminatod plate. Tsai^ developed a 
analytical method to predict the elastic constants of a lamina based on the 
properties of its constituent materials and extended the method to predict the 
elastic constants of the laminates. His work also described the application of 
classical lamination theory to laminated composite materials. 

The close form solutions to the governing differential equations for 
simply supported unbalanced cross ply and angle ply laminates under transverse 
load is presented by Whitney^. He concludes that the coupling effects becomes 
influential with increased degree of anisotropy and increased number of plies in 
the laminate. Whitny and Leissa^ presented a close form solution of simply 
supported laminated anisotropic rectangular plate using fourier series. 

Hottram and Selby^ have presented finite element fornvulation on the 
basis of assumed displacement. They have used 12— degree of freedom 
nor»onforming rectangular element. Jeyachandrabose and Kirkhop^ have presented 
finite element formulation using a triangular element and the stiffness matrix was 
derived explicitly by carrying out the basic matrix multiplication. 

There are several theories in the literature which account for 
transverse shear and normal stress, among them theories by Mindlin and Reissner 
are worth mentioning. Whitney and Paganc^^ extended the YNS theory to laminates 
coisisting of arbitrary number of bonded anisotropic layers. It was assumed that a 
normal to the mid plane before bending remains straight but not essentially normal 
to the mid plane, thus reclining use of a shear correction factor to transverse 
i^ear stiffness, the prediction of shear correction factor is problem dependent 
i.e. its value depeids on cmstituents ply properties, ply lay up, fiber orientation 
and the particular application. 

Chou and Carleone^^ presented a method for laminate analysis in 
i#iich the in— plane displacements are piece— wise continuous, the shear stress is* 
ocmstant tKJt cx>ntimiou5 across the laminate thickness and shear strain is 
disoontinumjs. Hintcn ja^escnted a finite element formulation which is based on 
assumed displacement model. He used 8— noded parabolic isoparametric plate 
bomiing' el«fient ' with tl-r'ee nodal degree of freedom. Pa.'nda a^nd Hatarajan^^ have 
used smer parametric 8— noded quadratic plate elem'^it with" 5— nodal' 'degree of 
freedcsm. R«KMy^^ used a C^ penalty plate 'bonding element for ' the solution of 
governing equations of' YNS theory. He solved a variety of ' problem's to 'Show the 



effect of material properties and the parametric effect of plate aspect ratio, side 
to thickness ratio, lamination scheme, lamination angle and number of layers on 
deflectim, stresses and variation frequencies of composite plate, the in— plane 
and transverse displacements and mid plane rotations were independently 
interpolated and coupled via the shear energy terms which behave as penalty terms 
that restricts the shear rotations to zero the thin plate limits. A finite element 
formulation based on the shear deformation theory was given by Lakshminarayan 
and Hurthy^^ usinof three— noded triangular element with 15— nodal dBgrees of 
freedom whioh includes displaooments and their derivatives. Engfolom and Ooho*^^ 
used equilibrium equations for a lamina and assumed parabolic variation of 
in— plane displacements and constant transverse displacement across the laminate 
thickness, to arrive at a more accurate prediction of transverse shear stresses. 
Kant and Sahani^^ used 9— node Lagrangian / heterosis element having five degrees 

JO 

of freedom per node. Akin and Kwon presented a mixed finite element formulation. 
They have used three moments and a transverse displacement as nodal degree of 
freedom. Iyengar and Umaretiya^® carried out deflection analysis of hybrid 
laminated rectangular and skew composite plates. They used the Galerkin approach 
to obtain results for kevlar / epoxy and boron / epoxy laminates. They concluded 
that for a given aspect ratio an increase in the skew angle increased the rigidity 
of the skew plate. They also concluded that for a specified deflection the hybrid 
laminates are lighter. 

20 21 

Pagano ' presented three— dimensional elasticity solution for 
composite laminates subjected to sinusoidal loads with simply supported bourwlary 
ccwiditions fcM' cylindrical and bi-directional bending problem. He concluded that 
classical lamination theory gave a poor prediction of laminate responses 
particajlarly at low span to depth ratio but converges to the exact soluticsi as the 
aspect ratio increases, the solution for the stress field converges mare rapidly 
than that for tha deflect ion. Later, Pagano and Wang^* further extended this work 
to include uiiformly distributed load and point load in the cylinck'ical bending 
IM'oblem. A hi#ier thecMry i^tfiich not only accoints for the transverse shear 

cteformaticns but also satisfied the zero shear stress conditicn at Uie top and 
bottom botr>ding planes was suggested by Reddy®® . Phan and Reddy®* proposed a 
finite element fcr'traalaticn in i^ich the in— plane displacemraits were expanded as 
cubic function of thickness coordinate and transverse displacement was kept 
ocMTstant. the contribution of transverse normal stress / strain was neglected. A 


mixed shear finite element based on higher theory was proposed by Putcha and 
Reddy** . They used eleven degrees of freedom per node , three displacements, two 
rotations and six moment resultants. Reissncr developed a theory by assuming 
the cubic variation of in— plane displacements and parabolic variation of 
transverse displacements across the plate thickness but neglected the terms 
contributing to the in— plane modes of deformations. He showed that the theory 
gave very accurate results when compared with the elasticity solution for the 
pure bending of an isotropic plate with a circular hole. Method of initial function 
has been used by Iyengar and Pandya ' to formulate the general problem of 
analyzing orthotropic rectangular thick plate as well as composite laminated 
plates. A differential equation of sixth order has been derived and a solution of 
the Navier type obtained to a particular case of uniformly distributed load on a 
siinoly supported square plate for various values of plate thickness and elastic 
constants. Pandya and Kant ' have formulated theories for generally orthotropic 
plates, which has been extended to laminated plates. Later Pandya and Kant have 
tmxlified the displacement model to include the in— plane degrees of freedom 
defined at the mid surface in the in— plane displacement fields. The vanishing of 
transverse shear ts^ms on the top and bottom bounding planes has been 
inccrporated in this formulation. 

1.3 Objectives and scx)pe of present work : 

Thou^ a multitude of papers have been published relating to higher 
crxier tl'reories , v&ry few of them incorporate finite element adaption . Also , the 
adaptability of a formulation for laminated composites , to the analysis of 
sarwlwich plates has been rarely tested. The responses predicted by higher cxxia' 
thew'ies require to be compared in order to establish their relative predictive 
capability. The theoretical assesment of responses have to be tested and cotwared 
with experimental results . 

The present work asseses the reliability and accuracy of suc^ a 
fonmjlation fcx' bending analysis of laminated plates under different physical 
coraiitions and using various material parameters. A hi^ner order theory frnaposed 
by Tarun Kant®* forms the basis for tha presant study . Tha FEM adaption of this 
fcxtnulation has been carried out by using C° , isoparametric, Lagrangian quadratic 
elenents . Transverse ^lear stresses have been evaluated using the constitutive 
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relations as well as the equilibrium equations . A large number of experiments 
were carried out to establish the reliability of this theoretical formulation and 
its FEH adaption. The results have been compared with those from the First Order 
Shear Deformation Theory and CPT in order to determine the domain of their 
applicability. A FORTRAN 77 program has been developed for the computer 
implementation of this formulation on the HP 9000 series 850 super mini computer 
systems in I I T Kanpar . 



CHAPTER 2 


THEORETICAL FORT-M-ATION 


2-1 Introductory Remarks : 

A Theoretical formulation based on the Higher Order Shear 
deformation Theory for the bending analysis of laminated composite plates is 
presented in this chapter . This theory accounts for realistic parabolic variation 
of transverse shear stress through the thickness . But it neglects transverse 
normal stress /strain effects . The governing equations , including expressions for 
interlaminar shear stresses are developed in terms of midplane displacements . 
Strain— displacement equations , constitutive relations and equilibrium equations 
are employed to obtain stress and strain values . 

2-2 Assumptions : 

( a ) Displacements and strains are small and in the linear elastic regeon . 

( fa ) The transverse normal stress and strain are neglegible . 

( c ) The points on the normal to the reference surface before deformation 

need not remain on a straight line after deformation . 

The assumption ( c ) is the basis for any higher order theory . It 
ifwlies a ncnlinear variation of the inplane displacements through the plate 
thickness . 


2-3 Definition of the Displacement Field : 


Taylor's series expansion in terms of thickness co-ordinate 'Z' is 
used for expressing the inplane displacements u ( x,y,z ) and v ( x,y,z ) at any point 
( x,y,z ) i^ile the transverse displacement w is kept constant through the thickness 
of the plate. In this way a 3— D elastic problem has been approximated as a 2— D 
plate px>blem . The displacement field is defined as follows : 


u ( x,y^ ) = u ( x,y/3 ) + Z 


V ( x,y,z ) = V ( x,y/3 ) + Z 


w ( XAJiZ ) — w ( xaj/3 ) 



(24) 


7 



where x,a,z represent the co-ordinate axes of the plate, ez = |^ = 0, 
that is this theory neglects the transverse normal stress /strain effects . The 
expressions for the inplane displacements u and v imply a nonlinear variation of 
these through the plate thickness thus the warping of the transverse 
cross-section is automatically incorporated . 

The various terms in the expressions 2.1 can be separated into two 
grcMJS representing membrane and flexural behaviour . 


membrane part flexure part 


u 

- Uo 

+ Z*u* 

+ 

Z ©X 

•f 

z" 

< 

V 

= Vo 


+ 

Z ©a 

+ 

z" 


w 

= 



Wo 





(2 -23) 


where Uq , vq are the inplane displacements and Wq is the transverse 
displacement of any point ( x,y ) on the mid— plane . and 0x represent the 
rotaticns of the normals to the mid— plane about Y and X axes respectively . The 
parameters ( u*, v* , 0* , 0* ) are the corresponding higher order terms in the 
Taylor's series expression and are also defined at mid— plane . Hence 


uo = u(x,y^);0x = [ ]□ ^ = fj 

Vo = V ( x,y,0 ) ; = [ §5 ]q = li 




ifu 

az" 
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lo 


Wo = w ( x,y^ ) 

The geometry of the plate along with the positive set of co— ordinate 
axes and the physical mid— plane displacement terms are shown in fig 2.1 . 

The meii±»rane parts can be dropped out from equation < 2.2a ) for a • 
s^jwnetric laminate subjected to transverse loading . This is because for a 
sv#Tm«tric laminate subjected to transverse loads , only the midplane of the plate 
r^a^senls a neutral str'face axvi there can not be any strain there . Hence for a 
s!#mnetric laminate , the sinwjlified expressims for inplane and transverse 




{x,y,z )- Laminate reference axes 


, FIG. 2-1. LAMINATE GEOMETRY WITH POSITIVE SET OF 

LAMINA/ LAMINATE REFERENCE AXES , DISPLACEMENT 
COMPONENTS AND FIBRE ORIENTATION. 



displacements become : 


U = z ©x + 7 ? e* 
V = z ey + z® e* 

W = Wo 


(2 -Zb) 


2-4 Strain— Displacement Relationships : 

The displacement field given by equations (2.2b) can be used to obtain 
strain at any point within a laminate . The theory of elasticity provides the 
definition of strain i4iich can be used to arrive at the linear strain— displacement 
relationship as follows . 


For uisymmetric laminates : 


£-x = ^ = €rXQ + Z Kx Hh Z* €xq Hh Z*^ Kx 

~ iy ^^0 + Z Ky + Z^ ejp -f Z® Ky 

'Txy = |y + = Cxvo + Z Kxv + Z* fcJv + Z® K^y 

'Tyi = + Z0a + Z® 

7 xz = ^ ^ + Z 0 y + Z® 

^i4iere 


C Kx # Ky j Kxy 


r ^0x 30x 

L 3x ' 8y ' 3y 




K 


xy 



ay ’ 8y ax _ 




(2‘3a) 


(2-3b) 


(2.3b) 




au* av* au* ^ 


( (^x . ) = [ 2 uj , 2 vj ] 


Strains for any layer 


shear strains and can be written io * 

[ei.] =[6o] +Z[K] +Z^[6;] ^ 

[ es ] = [ ^ ] + Z [ ] + Z* [ ^* ] 


Qan be grouped as boundary and transverse 
compact way as follows . 


where 



€x 


1 

X 

o 


Kx 

[^3 = 

ey 

7xy 

; [ £-0 3 = 

^yo 

T'xwo 

;[K] = 

Ky 

Kxy 


[‘;] = 


yo 


*ao 


;[k‘] = 


.* 

X 

* 

y 

.» 

xy 


and 


c £s 3 = [ ^ 3 + Z [ ^ ] + z"" [ ^* ] 


[<^3 


#y. 




th 




Similar c^E»rations for a synmetric case yi® 

[ efc 3 = Z [ K 3 + Z® [ K* ] 
and 

[ es 3 = [ ^ 3 + Z“ [ ^* ] 


(2 -48) 

(Z.4b) 


(2-58) 


(2.5b) 



2'5 Stress— Sti^ain Relations for an Orthotropic Lamina : 

For the sake of mathematical modelling the basic ply or lamina is 
considered to be a 2 —Dimensional , homogenous and orthotropic layer having 2 
principle material directions , one being parallel to fiber direction and other 
perpendicular to it . A laminate is formed by bonding laminae one over the other 
and is considered heterogenous through the thickness and generally anisotropic . 

Stress —strain relationships are obtained by taking into account the 
assumptions that w is constant throughout the thickness and that the transverse 
normal stress Oa = 0 . The generalized Hook's law for a homogenous orthotropic 
material can be written as follows . 

CTj = C;j fcj for i,3 = i , 2 , 3 , 4 , 5 , 6 

where Oj is the stress vector , Q; is the stiffness matrix and £; represents the 
engineering strain vector referred to the principal material axes 1,2,3 . 

Thus for our simplified case , the stress —strain relationship for an 
ortholropic lamina referred to the principal material axes can be written as 
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In a concise fcH'm it can be i^ittai as 

0-' = C t' (2-6) 

stress— strain relatiraiship can be referred with respect to the 


laminate reference system of axes ( x,y,z ) as follows : 

If a be the angle which the principal material axis i makes with the 
X-axis and principal axis 3 is coincident with the Z-axiSi then the direction 
cosines become : 



The transformation matrix T is given as 


Cos^a Sin^a ZSina Cosa 0 0 


Sin^oc Cos^a —ZSina Cos<x 0 0 


T = 


— Sina Cosa 


Sina Cosa Cos^a — Sin^a 


(2-7) 


0 0 0 Cosa —Sina 


0 0 0 Sina Cosa 


Transformation from one set of axes to the other is achieved using 

the relation 


Q = t"" C [t'1 


( 2 - 8 ) 


Her»ce the stress— strain relationships with respect to the laminate 



reference system of axes can be written as 
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Where coefficients of the Q matrix are defined as 


Qii = Cii Cos^a + ( 2 Cia + 4 Coo J Sin^a Cos*a + Caa SinV 

Qja ■■ ( Sin oc y- Cos ^ ) y* { ^ ^ Sin Cos a 

(iid = ( CjjL — Ci2 — 2 C33 ) Cosset Since + ( C12 ~ ^22 "f ^ ^33 ^ Sin^o: Cos<x 

Q22 ” ^11 Sin^cc + ( 2 Ci2 "f ^ ^33 ) Sin*<x Cos^oc + C22 Cos <x 

Q23 ^ i ^11 ^12 "“2 ^aa ) Sin^oe Cos<x -f ( Cj^a ^22 'f 2 Coa ) Cos cc Sinoc 

^33 "* i ^ ^12 ^ ^22 ‘ ” ^ ^33 ^ Sin oe Cos oe - 4 “ C33 ( Cos oc *4* Sin ot y 

Q 44 = C 44 Cos^a + Css Sin^a 

Q45 = ( Css - C44 ) Sina Cosa 

Qss = C44 Sin*a + Css Cos^a 

and Q|j •» , for i , j =» 1 to 5 . 


(2 * 10 ) 


2-6 Landnate ConstiUAive Relations : 

Tf» set of ecKJations expressing the flexure and the shear stress 
relaticr»sSiips for the differential element of a laminate as a fuBTctitn of mi<^lane 
curvature and shear rotations respectively , are knoi^r as the laminate 
constitutive relaticrts . 

The total iKJtential enerou IT of the plate with volumie V and surface 


area A can be written as 


ex' 


TT = U - W 

TT = I I&7 SL dV - E dA 
V A 


( 2 - 11 ) 


where U = strain energy of the plate 

W = work done by the externally applied loads 
£ = generalized displacement vector defined at the mid— plane 
E = Vector of force intensities corresponding to £ . 

Using the expressions for strains ( 2.3a ) and carrying out explicit 
integration tN'ouoSi the thickness of the plate we can get 

TT = I 2. dA - F dA 
A A 


i^ere 8 = , 6* / 8* ) 

For unsymmetric laminates e and a are defined as 


or 


* - [ ^0 ' ] 


arxi 


^ j^hJxjNy>Nxu»Njj^y^jjyJ^X»My^Xy,M^,My,Mj^y,Sx/Sy,TxiTy,S^,Syj 


(2.12a) 


(2.12b) 


or 


cr = L N , N . H . M , S . T . S J 



Nx 


n:' 

[ N ] = 

Ny 

.[n-] = 

n: 


Nxy 


< 


( 2 - 120 ) 



[M] = 


- 


“ 

Mx 



My 

1 1 

3C 

m 

V — 1 

It 

My* 

Mxy 


Mxy 


[S ] 


Sx 

Sy 


;[s‘] 


X 


;[T] 


Tx 

Ty 


where the components of the stress resultant vector cr are defined by 
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Similarly for symmetric laminates ^ 
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Stresses in equation ( 2.13 ) are now written in terms of the 
generalised strain components given by equations ( 2.12b ) and ( 2.14a ) by putting 
equation ( 2.4 ) or ( 2.5 ) , as the case may be, in the constitutive relations ( 2.9 ) . 

Thus the stresses in the i'th layer are given by the following 

relaticns : 

For unsymmetric laminates : 
u 

ay = [q j^eo ] + [q ]'" z [K ] + [q ]'’ 2 * [e* ] + [q ]V [k* ] (2-15a) 
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and for symmetric laminate 
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laminate constitutive relationship is obtained by substituting theese expressions 
fcM” stresses in equations ( 2.13 ) as follows 
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or of = D 6 (2'i7b) 

maiAw'ane rigidity matrix 


[ 
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membrane -flexure coupling rigidity matrix 
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flexure rigidity matrix 
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shear rigidity matrix 
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Similarly for symmetric laminates 
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The displacement model incorpcr'ating the above shear rigidity matrix 
is denoted by DM1 

The hi#ier order shear defcH'mation theory leads to non vanishing 
transva'se sdiear stress c»i the tc«D and bottom sur'face of the plate . To satisfy 
the condition of zero shear stress at top and bottom bounding planes ie. Ixz and 
equal to zero at z = dc ^ » the shear rigidity matrix is modified to incorporate 
these conditions as follows . 


= 0 are substituted in equations 


The conditions r7x*=7y*] 

( 2.5b ) to obtain 

_ if ( 2 - 19 ) 

These conditions are substituted in the shear rigidity matrix Dj such 
that the symmetric nature of Ds matrix is retained . 


Hence the modified shear rigidity matrix becomes : 
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The displacement model incoporating this modified shear rigidity 
matrix is denoted by DM2 


2-7 Expression for the Interlaminar Stress : 


The present theory accounts for the realistic parabolic variation of 
transverse shear stress through the plate thickness . But these strains when 
calculated by using the constitutive relations given by equations ( 2.4 ) and ( 2.5 ) 
violates the continuity requirement at the interface of two layers having 
non— identical Q matrices . Hence in this section interlaminar stresses are 
calculated by using ecRjilitr'iun equation of elasticity . 


For eat^ layer the equilibriixn equations can be written as 
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which can be rewritten as 
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The expressions for the in— plane stresses given by equations ( 2.5a ) 
and { 2.5b ) are used in the above equations ( 2. 2 la ) and ( 2.21b ) and integraticm 
is performed through the thickness of the layer . 

This yields for symmetric laminates : 
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CHAPTERS 


FINTIE ELEhENT MSCRETISATION 

3.1 PreliiBinarii Remarks : 

The dwnain is descretised into a number of finite elements ard a set 
of fuTctions is chosen to describe the state of displacement over each element in 
terms of its nodal values. The displacement function chosen enstres cc»npletr«ss 
within the element and cc»npatability across the inter— element bouTdaries. The 
total potential energy of the system ^hich is a fmction of assuned displacen»nt 
is calculated and tM^inciple of minimum potential energy is applied to gel the 
solutiiMTS of the fcrdanental equations of the higher order shear deformatimi 
theories f(»' laminated anisotropic plate, developed in chapter 2. Finite Element 
fcyrnwlation is carried exit for symmetric laminates only . Boundary Condiiicns used 
in this analysis is in fi^-ne 5.17 

32 Definition of displacemerA fuiction : 

The total solution domain is descretised into 'hEL’ siidomain 
(elements) such that 

NEL 

IT i6> = <6) C3.i> 


i^^rere n is tl-« potential of the system and If is the pot&itial of 
th® element . S is the vector of the urdenown displacement variable .The potential of 
an element can be i<«'itlen in terms of intemeal strain energy U® and the 
external mcxSc dona 14® as follows •' 

IT® (S) = U® - i32> 

For the displaaement model given by equation <2 2b) , the vectcB' of 
urtoxiwn displacement variable can be i«fl-'ilten as 

S - ( Wo , ©X ' ^ f <3-3> 

The same interpolation furKjtion is used to define all the components of 
gemralised displacement vector 6 such that 



in which hPE is the nurtoer of node element , Nj is the i^iape 
fuTctioo asscjciated with node i and Sj is the vKJtor of the generalise displacement 
associated with node i. In the present solution , the nine-noded isoparameta^ic 
Lagrangian quardilateral plate element is ccsTsidered whitdi has a paratolic 
variaticMn of the displaconent over the element. The shape functions are 
for ocM'Tier nodes i =i , 2 , 3 , 4 

^i = K ^ ^i ) ( ^ ^ ^i ) ( ^ ^i + ^ ^i “ ^ ) C3.5 a) 

for midside noctes i =5 , 6 , 7 , 8 

N® = |e €i ( 1+ € fi ) ( i- ) + I f? J?i ( i + T? f?i ) ( 1- I C35 b) 

for central ncKie 

Ng = { i- ) ( 1- ) (3.5 c) 


For iscK^arametric element the same ^lape function is used to 
describe the variation of unicnowi variable and geometric shape of elemait. f-tex>e 
to define tlie oxirdinate at any point on the reference xynalane of the plate we 






33 Elenent ^ffness Matrix : 

By expressing exterrtal strain enersu in terms of the unknown nodal 
displacanent we can obtain the stiffness matrix corresonding to the assjtned 
deformation state of an element .The strain en«^ is contirbuted n^iirane , 
flexure ar«i shear strain in the formulation -for uTsi^ranetric lam^te i^Tile in the 
case of ^jmmetric laminate manbrane , cxxitribution is a±>sa>t. 
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(3.8 a) 


(3.8 b) 


(3.8 c) 


(3.8 d) 


(3.8 e) 


Now using equations (3.4) and (3.7) , the generalised strain vectcB^ can 
be expressed in terms of mdal displacement as follows : 
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The internal strain energy of an element due to bending and shear 
action is obtained by integrating the porducts of flexural stress resultants with 
bending curvature and shear stress resultant with shear strains over the area of 
the element i.e. 
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The stress resultant can be substituted in terms of strains and 
curvature by using laminate constitutive relations given by equations (2.18). Hence 
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Substituting curvature and shear strain vector by the nodal 
displacement vector , using eqs. (3.9 a) and (3.9 b) we get 


U® = ||[ d’^(Bj Db Bt,) d + d'^(Bs Dg Bg) d ] dA 


= I d^|c| d 


< 3 . 12 ) 


where ihe stiffness matrix of the element is defined as 

[ K f = |[ (bJ Db Bt,) + (Bb Ds Bb) ] dA 


(3.13) 


3.4 Load Vector : 

Two types of ext^^al loading has been considered 
a) CcTK^eitrated rKKial load (Pc) , acting in the direction of the corresporKjing 
trodal degree of freedtwi 
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b) Uniformly distributed load (Pq) acting over the tcv face of the element in the Z 
directiwi. 

External woH< done by these fcw'ces on the element is given by 

<3.i4) 




Pc + 




Po dA 


= P® 


t^iere P is the element load vector 


Po 



_T 




<3.i5) 


where Pq^ is the vector of generalised concentrated load at node i in the 
direction of ccrresponding nodal degree of freeckwn. 


Po = (pIi. 

_T I 

where Poj: = ( Po 


^02' ^Oi' 

. 0 , 0 . 0 , 0 ) 



N = Matrix formed by ajatwxjpr lately placing shape function values. 


3.5 Derivation of the EquilibriLiR Equations : 


By eoB. <3.i) and <3Z) , total potential energy of the system can be 

exi»'essed as 


IT = 



c«i 


u® -w® 


By usir® eq. (3.12) ard (3.i5) , H can be i^itten as 
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d P 
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<3.i6) 


<3.i7) 


►Wch can te solved for the generalised displacenent 'd' 
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3£ Evaluaikm of the ^iffoess matrioes and Load Vectors : 


Gauss— Quartibature rule is used to evaluate all the integrals. The 
exact integratim of the integrals in the stifness matrices requires 3x3 Gauss 
point integration rule. This may lead to oversliff results at the thin plate limit, 
as the shear energy terms in the elastic stiffness matrix [K acts as p&^alty 
fuTKiticMi i<#uch cwTstrains the ^vz ^ zero as span to thidcness ratio is 

increased. To redi-Kse this effect a recfcjced e^z inte^ation rule is used for 
evaluating shear ene^ terms and 3x3 integration rule is adapted for ntx>-^iear 
terms. Hence el«nent of stiffness matrix is given by 
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and el^nent force vectcw' is given by 
1 1 


> P NG NG 

P® = Po + N Po df dT7 = Pc + ^2 N Po Wf Wg ! J 

JiJl ^ 


(320) 


where Wp and Wg are the weighting coefficients. 1 J 1 is the determirwnt of the 
Jacobian matrix relating the shape functicxn derivative in natural coordinate € and 
f? with those in X and Y coordinates. 



CHAPTER 4 


EXEPEJRIMEKrAL D'lVESTIGATION 

4.1 Material Qiaraclerization : 

Materials proposed for the use as facing or core of a sandwich 
construction were Alclad, bi-directional glass fabric, random glass fabric, cork, 
particle board, aluminium honeycomb, hard board and ordinary foam. Depending on 
their tensile modulus and shear modulus a final list of materials to be used for 
facing and core was prepared. 

4.1.1 Tests for material charecterization ; 

A. Tensile test; 

.Performed on the MTS 810.22 computer interfaced closed - loop 
servo - hydraulic material testing system for the facing 
materials 

.Performed on the INSTRON universal testing machine model 1195, 
for the core materials. 

B. Torsion test : 

. Performed on the INSTRON universal testing machine for the 
core materials only. 


42 Specimen Fabrication : 

42.1 Sandwich plate i^iecimen : 


Facing material 
Core material 


ALCLAD and bi - directional glass fabric 
A1 honeycomb, cork and particle board 
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Dimensions adopted 


Facing 

Core 

Length 

Width 

Facing 

thickness 

Core 

thickness 

ALCLAD 

Cork 

200 mm 

35 mm 

i.25 mm 

12 mm 

ALCLAD 

Part, board 

* $ 

# * 


* » 

Glass fab. 

Cork 

/ i 

t t 

i f 

f / 

Glass fab. 

Part, board 

$ > 

> > 

i t 

1 1 


Two rectangluar roselts were mounted on each specimen : one at the 
center of the bottom facing and the other at a diagonal location of 90 mm from 
the center in the two perpendicular directions. 

4.3 Experimental Set-Up : 

A fixture was designed and fabricated to simulate simply supported 
boundary condition on all the edges for the plate flexure test. It can accomodate 
plates up to 340 mm x 340 mm and 32 mm thickness. 

For applying uniformly distributed and concentrated load, two 
separate fixtures were designed and fabricated. All these fixtures were designed 
to fit into the grips of the MTS 810 servo - hydraulic material testing system. 

4.4 Test Procedure : 

Each specimen was secured properly in the test fixture and the 
readings of the various strain gauges were recorded with the help of a Multipoint 
measuring sustem. Care was taken that the load not exceed 30 % of the ultimate 
load for the facing to confine our investigations in the linear elastic range. The 
transverse deflectiwi was measured frcwn the movement of the cross head, and the 
strains were measured simultaneously with the help of a Multipoint measuring 
instrument. Tt® load was read directly from the digital display. 























CHAPTER 5 


RESULTS AND DISCUSSIONS 


5.i Prelimnary Remaks : 

The results obtained by using the present higher order formulation 
are presented and compared with the results obtained from various other theories 
available in the litrerature. The reliability , accuracy and applicability of the - 
present formulation and program for laminates with arbitrary geometry , loading, 
boundary conditions are examined through comparing its results with that obtained 
from Classical Plate Theory ( CPT ) , First Order Shear Deformation Theory 
< FSDT ) , other Higher Order Shear Deformation Theory ( HBDT ) and exact Three 
Dimensional Elasicity Solutions. 


Gauss quardrature scheme is employed to evaluate the element 
properties in the displacement based finite element formulation. The stresses and 
the stress resultants are computed at the gauss point nearest to the node 
concerned. The displacements are evaluaed at the node itself. Unless specified 
otherwise, 2x2 Gauss quadrature rule is applied for the shear terms and 3x3 
Gauss quadrature rule is applied for the flexure terms .The transverse shear terms 
are evaluated by using both equilibrium equations and constitutive relations. The 
superscripts 'ee' and 'cr' are used to represent the value of transverse shear 
stresses using equalibrium equations and constitutive relations respectively. 


The stresses , stress resultants and deflections are presented in 
non - dimensionlised forms as described below : 

Non-dimensional displacements (w,u,v) =KiX Actual displacement ( 

w , u ,v ) 


Non-dimensional inplane stressesCax , cry , Cxu) 


Non-dimensional transverse stresses ( ctxz , Oyi ) 


Non-dimensional bending resultants <mx,mu,mxw I 


Non-dimensional strains ( £xx > ^ 


Where 


Pn a^ 


K, = 


Poa' 


Ka X Actual inplane 
stresses ( Px < » o'xa > 

Ka X Actual transverse 
stresses ( Pxz • Pyx 1 
K 4 X Actual bending 
resusltantsC mx , my , mxu ) 
K 5 X Actual strains <£xx > 

£yV( } 
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Ks 


Poa 


K 5 = 


E_hi 

Po a 



( 3.1 ) 


Where ' h ' denotes thickness of the plate , ' a ' denotes a side of 
the plate and ' Pq ' denotes intensity of uniformly distributed load. 

The multipliers defined by equation ( 5.1 ) are modified for sinusoidal 
load by replacing Pq by Pc . 


Unless otherwise specified the location of various quantities 
presented are as follows 


u 

V 

Wq / mx I fTiy 
fBxy / CTxy 

0"x i CTy 

5x2 


Mid - edge of Y - axis at Z = 0.5 h 
Mid -- edge of X - axis at Z = 0.5 h 
Center of the plate at Z = 0.5 h 
Cornewr of the plate at Z = 0.5 h 
Center of the plate at Z = 0.5 h 
Mid - edge of Y - axis at Z = 0.0 
Mid - edge of X - axis at Z = 0.0 


5J2 Presentation of the Result and Discussions 

The results presented are classified into the following groups : 
a> Isotropic plates 

b> Orthotropic plates 

0 Multilayered sandwich plates 

d> Exeperimental results 


5JZ.i Isotrcwic plates 


A simply supported square isotropic plate subiected to uniformly 


distributed transverse load is considered with the following properties - 

3 

Ei = E2 = E ; ^12 = ^21 = 0.3 ; 0]__ = i 


12 ( 1 




0x2 ~ ^23 = G 


'23 


^13 


2 < i 4- 1 ; I 


a 


a * 10 ; 


Pq =, IfM 


a / h = 4 # 5 , 10 and 100 . 



Results of convergence of maximum deflection , stress resultants and 
stresses for the plate with side to thickness ratio equal to iO are presented in 
Table 5.1 for the two displacement models DMi and DMa . Table 5.2 contains the 
values of maximum deflection , stress resultants , and stresses for a/h = 5 and 
100 . Results available from other higher theories , Elasticity solutions and 
Classical plate theory ( CRT ) are also presented in these tables for comparision 
purposes . Parabolic distribution of transverse shear stress through the plate 
thickness is obtained by higher order shear deformation theory as depicted in 
fig. 5.1 . The variation of normal inplane stress through the plate thickness is 
presented in fig . 5.2 for simply supported and fixed boundary conditions. The 
effect of plate's side to thickness ratio on the warping of the transverse cross - 
section is shown in fig . 5.3a through 5.3b . 

The following observations are made from these results : 

• For the thick plate ( a/h = 5 ), the results obtained from different sources 
vary marginally but they all converge to the classical thin plate solutions in the 
thin plate domain ( a/h = 100 ) . 

• The warping of the transverse cross - section is more pronounced in the 
case of thick plate as is seen from figs. 5.3a through 5.3b . 

• From the Tables 5.1 , it is observed that the values of the transverse 
shear stress obtained from constitutive relationship are closer to the Elasticity 
solution as compared to that ctotained by equilibrium equations . 

• 3 x3 mesh in quater plate is sufficient to predict deflections and in - 
plane stresses but , more refined mess is required for accurate prediction of 
stress resultants and transverse shear stresses . 

52JZ Orthotropic plates : 

Example 1 : 

A foiT’-ply square laminate < 0/90/90/0deg. ) with layers of equal 
thickness and subjected to sinsoidally distributed transverse load. 

The lamina propa'ties are assumied to be 
Ejl “ 25 X 10 psi i Ej “ 10 psi r ^12 ** ^13 0.5 x 10 pssi j 

G 23 =» 0.2 X 10 *psi ; i>i 2 <- “ 0-25 > 

Table 5.3 contains the value of non - dimensionalized deflecticrvs and 
stresses . The tlr^e dimwisional Elasticity solutions of Pagano and Hatfield for 



simply supported rectangular plates under sinusoidal loading are used to asses 
the improvement . Figure 5.4 depicts the effect of material anisotropy on 
transverse displacement . 


The various quantities have been nondimensionalized as follows 


w = 


w h'’ E, 


Po 


X 100 


; w = w < a/2 , b/2 ) ; 


= ( a/2 , b/2 , h/2 ) x-O-^ 

Po a“ 


; = CTy { a/2 , b/2 , h/6 ) x 


P 2 
0 a 


= CXxy (0,0, h/2 ) 


h . 


Po a 


Cyz = ayz ( a/2 ,0,0) p — ^ ; 


o'xz = (Jxz ( 0 , b/2 , 0 ) i 


From the table 5.3 and figure 5.4 , the following conclusions are 

drawn : 

• The exact stresses Hx , Cy and ctxu computed using the constitutive 
equations of the higher order theory , are greatly improved over the the results 
obtained using the first order and the classical plate theory . 

• From fig. 5.4 it can be concluded that the intensity of shear deformation 
also depends on the material anisotropy of the layers . The maximum transverse 
displacement at center reduces considerably as the degree of orthotropy 
increases . The classical plate theory underpredicts the deflection even at lower 
ratio of moduli . THe disparity between the higher order and the frist order 
theory may be due to the higher order contributions of the present theory and 
shear correction factor whose value depends on the lamina properties and the 
lamination scheme . 

Example Z : 

A singly siworted suqare orthotropic plate subjected to tniform 
transverse pressLH'e is condsidered with the following properties : 

= 3 , 10 ,40 ; = ^ = 0.6 ; ^ = 0.5 ; 

Ey Ey Ey Ey 

Ey ; T*y « 7y* - Ixx - 0.25 ; Po -iOO ; « - 0° ; 
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The maximum transverse displacement for a/hs=3 ,10 and 40 have been 
shown in table 5.4 and the results available by first w'der shear defcM'mation 
theory , Reissner's theory and method of initial function have also been 
represented for compariscwi sake. The variation of in-plane displacement <u & v) 
has/e been depicted in figure 5.5 . The variation of in-plane normal stress and 
transverse shear stress has been shc»«i in fig. 5.6a and 5.00. The following 
observations are made from the results : 

• That the present formulation incorporates the warping of the transv®-^ 
cross-section is evident frtwn fig. 5.5 . 

• The maximum transverse displacement at centre reckices considerably as 
the degree of orlhotrc^jy increases as shorn in fig. 5.7 

• The maximten discreparxoy in di^lacemeit of about ZS*/. is observed with 
displacement model DM2 for a/h s= 3 and Ex/Ey = 10 . Otherwise , the preseot 
results differ OTtly marginally from the available solutions for various Ex/Ey and 
a/h ratios . 

5JZ3I Multilayer Soidwich Plates : 

A sirw>ly suwxjrted fi\«— layered square isotropic sanci^iic^ plate 
under urtifonti trarsverse fH^ssur^ is considered with the following material 
properties (as shon in Fig. 5.13) : 

Case — la 

Stiff Layers : E=1Q^ psi ; h=0.02 in. Cw^s; 0=3x10^ psi ; i'=0.0 ;c=0.4 in. 

Case — Ib 

Stiff Layers -• E—10^ psi ; i»«=0.3; h>»0.02 in. Cores: G^SxlO® psi ; i/«0.0 x5=0-4 in. 
Case — Ila 

Stiff Lay®^ San« as in case la 

Ceres Same as is used in the case la except that the the bending 

rigidities of the core are neglected . 

Case - Ilb 

Stiff Lay»^ Same as in esase Ib 

Cores Same as is used in the ca^ Ib exo^ that the the bending 

rigidities of the core are reglected . 
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Case — Ilia 

Stiff Latr^es Same as is used in tlie case la except that the shear rigidities of 
facings are neglected . 

Cores Same as is used in the case Ua 

Case - Illb 

Stiff La^^s :Same as is used in the case Ib except that tt« shear rigidities of 
facings are neglected . 

Cores Same as is used in the case lib 

NOTE Fot' the sake of conwarison: 

The results of FOST for case I are given fc»' tref . 323 

Stiff Lagars psi ; ii'iBO.3; h«»0.02 in. 

Core below the mid - plane : G»3xi0^ l»i ; t/ » 0.0 ; o » 0.4 in. 

CcM^ above the mid - plane : G=5xi0^ psi ; v = 013 ; c *= 0.4 in. 

The results of FOST ior case II are given fo' Cref . 323 

Stiff Layers Same as in case I 

Cores Same as in case I except that the beiding rigidities of the cores 

are neglected 

The reeults of FOST, FEX** and Series Solution for case ni are given 

for tref. 323 

Stiff Layers Same as in case I exceat ^t the shear rigidities of facings are 
mglected . 

Ccw'es Same as in case I exc^t that the bending rigidities of the cores 

are neglect«i 
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The stress resultants and central maximum transverse displacenent for the three 
cases have been shoMn in table 5.5 .The follc»dng enervations are made 

• The effect of ccH^idering bending rigidity of core is only marginal as is 
seen clearly by comparing cases I and II. 

• The theories t<»rfTich do not cwisider shear rigidity of stiff layers of the 
order 3.8xl0*psi and bending rigidity of core of the order l.'^xlO^psi over 
estimates the deflection by a very large margine as conpared to the thecries 
which consider bending and shear rigidities of both facing arwi core. 

53 Experimenial Results : 

53J. Sandb>iiich paiels : 

The values of Yojtg's mtxiulus and shear mcxiulus obtained from 
tensile and torsion tests are presented in Tables 5.5 . The values of v adopeted 
for various materials are also ireiicated in the same table . Fixtures for holding 
plates arxi acplying ipifconly distributed load are shc»4n in fi^jres 5.14 to 5.16 

The values of maximum transs/erse deflectios ( measured at the center 
of the botton facing ) , maxiimjm normal irplane strains (also measured at 
center of the Ixjttom facing ) and maxiim^ inplane ^lear stresses ( measured at the 
ccM'ner of the bottom facing ) are Fx'esented in figures 5.8a to 5.11c . The 
theoretical valies of maxinttjn trarsverse displacements , maxiimm) inplane normal 
strains and maximLan inplane shear strains obtairpd by i^ng the present FEM 
fonmjlation and by Alibi's®’' formulation are also depiced in the same fisaJres . 

It is dtserved that experiir«ntal results for most of Wie specimens 
cwif arm with the thecM~atical fpedictions . Moreover , theses values lie close to 
the theoretical values for all ^jecimens at anall strains . The expe'in^ital values 
lie close to the predictions of the pre5®tt FB4 f ormilaticns as cewpared to that 
of Allen's . At higher str^s-strain levels ext^riraental values cteviates from the 
thecx^tical predictions . This may be due to tt® cause that at the higher stress 
and strains levels the materials may becom rexvlinear elastic or even plastic , 
violating the lineer elastic assimptions in «« abo\« mentioned theoratical 
fpmilatiors . 



Convergence of maximum deflection, stress resultant and stresses in a simply supported square isotropic plate 

under uniform transverse pressure (a/h=iO) 
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Maximum deflection, stress-resultant and stresses in a simply supported 
square isotropic plate under uniform transverse pressure 
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Nondimensionlized deflections amd stresses in a four- layer cross ply (0/90/90/0 deg) square 
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Table SA 

Maximum transverse displacement in a simply 
supported square orthotropic plate 
under uniform transverse pressure 

Ex / Ey 

Source a/h 3 10 40 

DM1 5 .3448E-0i .2198E-01 .i400E-01 

DM2 .3711E-01 .2529E-01 .180iE-01 

FOST .3445E-01 .2201E-0i .i422E-01 

Method of ini- 

-tial functions .3445E-01 .22i5E-Ql .i472E-0i 

Reissner's 

theory .3389E-01 .2172E-01 .1410E-01 

DMi 10 .2948E-0i .1585E-0i .6752E-02 

DM2 .3014E-01 .1672E-01 .6844E-02 

POST .2945E-01 .1584E-04 ,6760E-02 

Method of ini- 

-tial functions .2938e-0i .1582E-0i .6780E-02 

Reissner's 

theory J2933E-0i .1578E-01 .675Cffi:-02 
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Table 5S 

Maximum displacement and stress-resultant in a 
simply suppc»Hed five-layered square isotrcxoic sandwich plate 
under uniform transverse load 


Source 

Case 

Wo 

ff>x 



la 

.3480 

.02272 

01466 


Ib 

.3609 

.02333 

01496 

POST* 


.34^ 

04750 

.03103 

mi 

Ha 

.3392 

X32334 

01320 


Hb 

.3628 

02347 

01504 

POST* 


.3555 

.04820 

.03160 

I»11 

nia 

.5286 

02435 

OlKl 


nib 

.6702 

02613 

01998 

FCKT* 

1 

.6669 

0482 

03160 

FEM®« 


.8166 

0446 


^RES* 


8360 

— 



SOLUTION 






Case 1 : CcMTsidering bexling and shear rigidities 

for toth facings and ctw'es . 

Case II : Considering bending and bt^z’' rigidities 

for facings Ixit wily transverse shear 
rigidities for cores . 

Case HI : CwTSic^ing only beiding rigidities fix' 

facings and trarsverse stwar rigidities for 
wwres . 


• Fwr stawiard results please see note on page ^ 




Table 5 jB 

Estimated Value erf Yotsig's Hociulus 3hear Mcxlulus 
arwd Poissen's Ratio for different materials 


Material 

Youig's Modulus 

E (GPa) 

Shear Mcxkjlus 

G(MPa) 

Poisson's Ratio 

u 

Alcalad 

7225831 

2779i£58 

05 


14.82812 

K;40.9545 

QZl 

Cork 

0J3131528 

45714285 

055 






PBoard 


0.1571288 


59.518518 


0 .^ 























Fig. 5.2 Vonation of normal inplone stress through . j 

the ploie thickness for a simply supported d.nd rixecf 
isotropic plate under uniform tronsversa lood 




Fig. 5.3b variation of ini^iane displacements thruogh 

the plate thickness for a simply supported 
isotropic plat© under uniform pressure . . 
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Fig, 5.7 Effect of degree of orthotropy (Ex/Ey) on the 
maximum transverse displacement 
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Maximum in— plane shear strain at the bottom facing of 

a simply supported sandwich panel with Gloss fabric foclng 
and particle bound core under unifonmolly distHbuted load 




(b) STRESS-RESULTANT COMPONENTS ( FORCES ) 



(c) STRESS- RESULTANT COMPONENTS (COUPLES ) 


FIG. 5-12 POSITIVE SET OF LAMINA STRESSES AND LAMINATE 
STRESS RESULTANTS (FORCES AND COUPLES) 
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FIG. 5.13 FIVE LAYERED SANDHICH PLATE FOR SECTION 5.2. 3 











CHAPTER 6 


CONCLUSIONS 

6.1 toieral : 

A higher order shear deformation theory for bending analysis of 
symmetric laminate has been developed . Extensive experimentations on various 
sandwich beam and plate specimen were carried out . On comparing the results 
obtained from FEM analysis, experimentation and various other classical theories , 
first order shear deformation theories and higher order shear deformation 
theories , the following conclusions can be drawn : 

• The higher order shear deformation theory accounts for parabolic 
distribution of transverse shear stresses through the plate thickness and hence 
it does not require shear correction coefficients. The first order shear 
deformation theory use arbitrary shear correction coefficient) which is not 
justified as the value of this shear coefficient depends on lamina material 
properties and stacking sequences. 

t The classical laminate theory predicts correct responses only in thin plate 
domain (a/n > 20 >. Though the first order shear deformation theory gives 
satisfactory results for transverse displacement, but it fails to predict corhect 
value of stress for thick laminate. But the higher order shear deformation theory 
can be relied upon for deflection and stress calculation for thin as well as thick 
laminates. 

• The higher order shear deformation theory shows considerable warping of 
the transverse cross-section, especially in the case of thick laminate with weak 
core. This type of realistic behaviour cannot be predicted by first order or 
classical shear deformation theory. 

• The ordinary sandwich plate theories neglect the shear rigidity of the 
facings and the bending rigidity of the core. They over estimate the deflecticn. It 
is seen here that bending rigidity is not much influence on the result , but if 
shear rigidity of facing is added to the bending rigidity , then there is 
considerable change in the deflecticm. 

62 Reconmienrfations for fl^xjre moHc : 

The fM^sent investigation is amenable to further extension or 


improvement as suggested below : 

• The present formulation for the plate structure can be extended to deal 
with composite or sandwich shell structure. Also , it can be extended to stiffened 
plates and shells. 

• The laminate plate or sandwich structure can be subjected to temperature 
gradient across its thickness, as the outer exposed surface may be sub'^cted to 
the hot jet blast or aerodynamic heating. The present Theory can be revised to 
incorporate such effects. Also the experimental investigations can incorporate 
simultion of more real life settings and a data base can be developed. 
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